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ABSTRACT 



We study thermal emission from circumstellar structures heated by gamma-ray 
burst (GRB) radiation and ejecta and calculate its contribution to GRB optical and 
X-ray afterglows using the modified radiation hydro-code STELLA, ft is shown that 
thermal emission originating in heated dense shells around the GRB progenitor star 
can reproduce X-ray plateaus (like observed in GRB 050904, 070f fO) as well as de- 
viations from a power law fading observed in optical afterglows of some GRBs (e.g. 
020124, 030328, 030429X, 050904). Thermal radiation pressure in the heated circum- 
burst shell dominates the gas pressure, producing rapid expansion of matter similar to 
supenova-like explosions close to opacity or radiation flux density jumps in the circum- 
burst medium. This phenomenon can be responsible for so-called supernova bumps in 
optical afterglows of several GRBs. Such a 'quasi-supernova' suggests interpretation 
of the GRB-SN connection which does not directly involve the explosion of the GRB 
progenitor star. 
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1 INTRODUCTION 

Cosmological gamma-ray bursts (GRBs) are the most lu- 
minous events in the Universe, which are observed as short 
flashes of gamma-ray radiation (prompt emission) accom- 
panied by tra nsient afterglow on longer wavelengths (see 
iGehrels. Rami rez-Ruiz & Fox 2003 for a review and refer- 
ences) . 

After decades of intensive studies, it has been commonly 
accepted that GRB radiation has mainly a non-thermal na- 
ture. The principal prompt emission mechanisms include 
synchrotron radiation of electrons accelerated in relativistic 
shock waves or magnetohydrodynamic (MHD) instabilities, 
synchrotron self-Compton and inverse Compton scattering 
of thermal photon s (see e.g. lBeloborodovll2010l : lDaigne et al.l 
l201ll : IZhanj|201ll ). 

However, thermal emission can also emerge in GRB 
light curves for several reasons. It is commonly accepted 
that long GRBs are connected with explosive deaths of mas- 
sive stars, which can be surrounded by a dense stellar wind 
or even by shells of matter (up to several Mq) produced 
by pulsational instability of the progenitor. Even several 
per cent of the huge GRB energy ~ 10^" — 10^'' erg inter- 
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cepted by this matter can b e observed as thermal emission 
l|Blinnikov fc Postnov|[l998l ). 

The rmal emission si gnatu res have been found in GRB 
spectra. iRvde fc Pe'erl (|2009l ) calculated a photospheric 
component on top of the non-thermal power-law, which 
can describe several tens of observed GRB pro mpt emission 
spectra better than a pure broken power law IjBand et al.l 
1993) or a cut-off power law, i.e. the 'can onical' spectral 
functi ons. Several detecti ons of X-ray lines (lAntonelli et al.l 
I2000| : Reeves et al.l 20021') and blackbody c omponents (see 
e.g. ICampana et al.l |2006| , iPage et aU I2OIII ') were reported 
in X-ray afterglows suggesting that thermal effects may also 
contribute to the afterglow radiation. 

iBlinnikov. Kozvreva fc Panchenkd (|l999l ) dis- 

cussed how the power-law GRB spectra may be 
produced by the blackbody radiation with chang- 
ing temperature. Recently, studies by several au- 



thors (e.g. Lazzati. Morsonv fc Begelman 



Mizuta. Nagataki fc Aoil I2OI1T iNagakura et al.l 



200g, 



2011 



Suzuki fc Shigevamall2013l ) make the GRB thermal prompt 
emission more topical. Using axisymmetric relativistic 
hydrodynamical simulations, they argue that, breaking 
through the progenitor star envelope, the relativistic jet 
becomes coUimated and cut apart by tangential shock 
waves. Different parts of the jet start moving with different 
Lorentz-factors (an order of magnitude exceeding those 
predicted by the spherically symmetric fireball model) and, 
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when colliding with each other, they can produce thermal 
emission. This l ooks similar to the widel y accepted in ternal 
shock scenario (JRees fc MeszarosI Il994l : iPiranI Il999l ') , but 
the emerging thermal emission ensures greater efficiency 
than the synchrotron one. 

In addition, late-time afterglow light curves of GRBs 
often demonstrate bumps with colour evolution and op- 
tical and NIR spectra looking like type Ib/c supernova 
on top of the power-law fading. These supernova-like fea- 
tures have been detected photomet rically and spectroscopi- 
cally in m any GRB afterglows (see lWooslev fc BloomluOOq . 
ICano et al...201L for a review of GRB-SN physical and obser- 
vational connections), and are considered as the most com- 
pelling evidence for the presence of thermal plasma radiation 
in GRB afterglows. This GRB-SN connection is a corner- 
stone of the present concepts of the origin of GRB central 
engine. 

Therefore, thermal emission in GRB afterglows can ap- 
pear as irregular, 'sporadic' deviations from the 'canonical' 
power laws of GRB light curves and spectra F^, ~ t~°'v~^ . 

Thermal emission from large circumburst 
structures has also been studied. For ex- 



ample 



iBisnovatvi-Kogan fc TimokhinI (|l997h : 



iBarkov fc Bisnovatvi-KoganI (2005) studied the case 



where a dense molecular cloud around and near the GRB 
progenitor star was present. In their simulations, a cloud 
with number density n ~ 10^ cm~"^ and radius of a few 
parsec was heated up by an instant (delta-like) gamma-ray 
emission pulse via Compton scattering, and the generated 

thermal radiation was found t o last for a few years. 

In our previous paper (jBadiin. Blinnikov fc Postnovl 

l20ld . hereafter Paper I), we modelled thermal effects in 
GRB afterglows using the multigroup radiation-hydrocode 
STELLA. The code was originally dev eloped for detailed sim - 
ulations of supernova radiation (see iBlinnikov et al.lll99a ). 
and is suitable for calculation of thermal emission generation 
and transport in exploding and expanding media. In paper 
I the code was modified to include gamma-ray heating and 
non- stationary ionization of matter. 

In the present paper, the code is modified in order to 
more accurately handle discontinuities (in degree of ioniza- 
tion, opacity and emissivity, etc.) which emerge due to the 
propagation of gamma-ray radiation. The GRB jet inter- 
action with the shell is modelled using the so called 'ther- 
mal bomb' (in contrast to a heavy 'piston', or 'quasi-ejecta' 
which we used in Paper I). In this approach a certain amount 
of thermal energy is deposited into a certain amount of mass 
during a certain time interval. Different boundary conditions 
can be imposed to thermal radiation escaping from the shell. 

In Section 2 we describe the model. Section 3 presents 
results of calculations. In Section 4 we compare the results 
of our calculations with features of observed afterglows and 
discuss other possible astrophysical applications. Section 5 
summarizes the results obtained. 



the cold dense shell surrounding the GRB progenitor and to 
find possible observational signatures of this radiation. 

A massive GRB progenitor star can experience a pow- 
erful pulsational instability stage before the explosion and 
expell severa l M^ of matter. This scenario , for ex ample, was 
suggested bv lWooslev. Blinnikov fc Hcgcr (2007) to explain 
an extremely bright Type Iln supernova SN2006gy. Observa- 
tional signatures of such mass-loss events arou n d sup ernova 
progenitors were recently found bylOfek et al.l (|2013l ). 

Like in IWooslev et al.l ((20071), in our calculations we 
used the multigroup radiation hydrocode STELLA. The code 
is one-dimensional and non-relativistic, but radiative trans- 
fer equations include aberration, time delay and Doppler- 
shift effects to within the 0{v/c) accuracy. Therefore, to 
aviod relativistic motions, we are restricted by spherically 
symmetric shells with densities corresponding a total Thom- 
son optical depth of less than one and located far enough 
from the centre of explosion. Specifically, we consider thin 
shells at a distance of ~ 10^® cm with number densities of 
several 10^" cm" 



As IWooslev et al.l (|2007|) show, such and even denser 
'walls' around the GRB progenitor could arise if the central 
star had experienced several mass ejection episodes, and the 
ejecta interacted with each other. 



2.1 Initial parameters 

Let us consider a spherical shell with mass 5 Mq, radius 
R — 10^^ cm and thickness 5x10^'^ cm. The mean baryon 
number density in the shell is ribar ~ 10^^ cm"'', which 
is high enough to produce considerable thermal effects and 
keep the dynamics non-relativistic. The relatively small shell 
thickness allows a good spatial resolution. The elemental 
abundance was assumed to be the s ame as in the o u ter sh ells 
of the pre-supernova calculated bv lWooslev et al.l (|2007|). 

The shell is illuminated by gamma-ray radiation with 
a peak isotropic luminosity of 3x10^'' erg s~^ consisting of 
three FRED (Fast Rise - Exponential Decay) pulses, each 
with the characteristic duration 1.3 s, so the total energy of 
gamma-ray emission was about 4.5x10^'' erg. These lumi- 
nosity and the gamma-ray peak duration are quite typica l 
(see, e.g., statistical study in lBadiin. Beskin fc Grecdl2009l V 
The gamma-ray light curve was intentionally made 'spiky' 
to see the possible quick response of the surrounding shell to 
such rapidly altering 'perturbations' (no effect was actually 
found). 

The spectral energy distribution of the gamma-ray 
emission was taken in the form of a broken power law (the 
Band function, Band et al. 19_93,) with the characteristic en- 
ergy eo = 300 keV and low- and high-energy exponents 0.9 
and 2, respectively. The half opening angle of the gamma-ray 
emission cone was taken to be Sjet = 10° . To treat the emer- 
gent thermal emission, 120 frequency groups from 50000A 
in IR to 30 keV in X-rays were used. 



2 THE MODEL 

The primary goal of our calculations is to study ther- 
mal emission generated during the interaction of powerful 
gamma-ray fiux and following blast wave from a GRB with 



2.2 Abilities of the code 

For the present calculations several modifications were im- 
plemented to the STELLA code. To treat the interaction of 
gamma-rays with matter, a time-dependent system of ion- 
ization balance equations is solved. Photoionization, ioniza- 
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tion due to Compton scattering on bound electrons, coUi- 
sional ionization, radiative and dielectronic recombination 
are taken into account. These elementary processes deter- 
mine the effect of the gamma-ray emission on the state of 
matter, its opacity and emissivity. 

The opacity of the shell matter should not block the 
non-thermal synchrotron emission of relativistic shocks in 
the GRB ejecta (in order to observe the GRB phenomenon 
itself). This requirement is fulfilled because the full ioniza- 
tion of matter by gamma-rays strongly reduces its photo- 
electric opacity. A large value of the ionization parameter 
^ = 978nj^Q^ for photon energies above 10 keV and 1760 n^p^ 
for those above 1 keV (nio is the bound electron number den- 
sity in units 10^" cm~^) indicates that there is a lot of high 
energy photons per every bound electron. Our numerical 
non-stationary ionization calculations show that the matter 
becomes fully ionized in the first thousandths of a second 
after the gamma-ray forward front. Only a few hydrogen- 
like ions of heavy elements (e.g. FeXXVl) with the high- 
est recombination rates remain in non- negligible amounts of 
0.1 — 10 cm~^. Thus, both the prompt non-thermal GRB 
emission and the early X-ray afterglow will be visible until 
the recombination of matter occurs. 

The radiative heating of cold matter is produced by 
the gamma-ray radiation via photoionization and Compton 
scattering (with opacity coefficients depending on the non- 
stationarily calculated state). In addition to the heating ef- 
fect, the radiation drag force produced by gamma-radiation 
is taken into account. 

As the radiation front moves with the speed of light, 
the unhealed regions are not causally connected to heated 
ones; there is a sharp discontinuity in temperature, ioniza- 
tion, opacity and emissivity behind the gamma-ray front. 
No outward fluxes of heat, mass and radiation through this 
contact discontinuity are assumed as boundary conditions. 
In the modified code this is taken into account. When cal- 
culating the Eddington factors of the thermal emission (see 
[Blinnikov ct al. 1998), the heated and unhealed parts of the 
shell are treated independently. 

Due to high densities and radiation fluxes (both thermal 
and non-thermal), the radiation pressure can dominate the 
gas pressure by four orders of magnitude. This makes the 
problem very sensitive to the local opacity, direction and 
strength of radiation flows. Depending on whether there are 
highly opaque regions (where the radiation is strongly cou- 
pled with matter) or how transparent the inner boundary 
is (we mention here only the inner boundary, since the flux 
conditions on the outer and intermediate ones are set strictly 
as 'outwards only' and 'inwards only'), the radiation either 
flows away exhausting shell's thermal energy, or is accumu- 
lated in hot regions and drives their expansion. 

The latter, for example, takes place if the inner bound- 
ary for some reason does not transmit the radiation towards 
the centre (e.g. it is highly compressed by a shock and rather 
opaque, or scatters the radiation back, or due to relativis- 
tic aberration makes most of the emission flow outwards). 
Then the thermal radiation appears trapped between the 
'no-flux-towards-the-centre' inner and 'no flux outwards' in- 
termediate boundaries, accelerating the hot zones rapidly up 
to several 10* km s~^, which results in a 'quasi-supernova' 
event (see Section 4) . As matter is needed to be compressed 
by several orders of magnitude to compensate the radiation 



pressure, a radiation implosion easily develops around the 
opacity or radiation flux jumps. Note that the compression 
increases the opacity, providing a positive feedback on the 
implosion, until the radiative pressure growth coasts or some 
instabilities come into play smearing off the density peak. 

We place the inner boundary flux condition by means 
of the inner boundary Eddington factors. They are calcu- 
lated similarly to those at the outer bound ary in the origi- 
nal STELLA code {Hu = h^Jv in terms of iBlinnikov et al.l 
:199g). An arbitrary incoming boundary flux, e.g. a diluted 
blackbody background or a counter-jet emission, can also be 
taken into account. 

In addition to being heated by gamma-rays, the shell is 
heated due to partial thermalization of the GRB jet kinetic 
energy. Clearly, the one-dimensional non-relativistic code 
cannot model the relativistic jet dynamics self-consistently; 
therefore we model the jet impact by means of a 'thermal 
bomb'. In this approach, which is widely used in simulations 
of supernovae, the jet braking and its kinetic energy dissi- 
pation are described in terms of the thermal energy which is 
deposited into selected zones during the time interval mostly 
appropriate for the effecive dissipation process. 

This general approach is rather insensitive to details of 
the GRB flreball expansion or blast wave propagation. As- 
suming energy equipartition between the gamma-rays and 
ejecta, we have taken 4.5 x 10^^ erg as an estimate of the 
energy to be thermalized. In fact, results of calculations 
turned out to be sensitive to the total thermalized energy 
and weakly depend on deviations from the equipartition. 

In order to determine where to deposit the energy in the 
shell, the shock Lorentz factor should be known. The 'ther- 
mal bomb' domain should have mass exceeding Es\i{cTs\-i)~^ 
[the characteristic swe pt mass scale at which the relativis- 
tic self-similar regime (IBlandford fc McKeelll976r ) is estab- 
lished], but less than i?shC~^ (the characteristic scale of 
the non-relativistic Sedov- Taylor regime). The blast wave 
moving with the Lorentz-factor T^h impacts the shell in- 
ner boundary at the time (Jt^-sh ~ Rs\-ic\\{'icV'i^)~^ after 
gamma-rays. The duration of the heating process can be es- 
timated as the radial thickness of zones in which the energy 
is deposited divided by the speed of light. 

We assumed the blast wave Lorentz-factor Fgh ~ 30. 
This is suitable for the incoming shock with the initial 
flreball Lorentz-factor Fsh.i ~ 10^'^"'' which follows the 
Blandford-McKee solution with Fsh oc R~^'^ in a homo- 
geneous medium for an adiabatic shock (or Fsh oc i?""^ 
for a radiative shock) well before the impact. This means 
that several -Bsh(cFsh,i)"^ ~ 2.5 x lO"'' - IO'^'Mq and 
^ Es\iC~'^ « 0.25Mq (in order for the shock to be rela- 
tivistic) of circumburst matter should be swept up by that 
moment. This indeed could be the case if the shell itself were 
produced by a non-relativistic shock initiated by one non- 
relativistic ejecta mo ving through anothe r (like in a pair- 
instability supernova. IWooslev et al.ll2007l ). 

Thus, with Fsh ~ 30 we have placed the 'thermal bomb' 
within 20-Esh(cFah)~'^ ~ 0.05Mq of shell's matter from its in- 
ner boundary. This also corresponds to the first Lagrangean 
zone. The 'thermal bomb' is 'ignited' 5t~f-sh « 200 seconds 
after the radiative heating begins in the first zone and is ac- 
tive for 17 s (which also coincides with the prompt emission 
duration). Hence, there are two kinds of heating processes 
involved: the radiative heating which has the 17 s transient 
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profile and moves outwards from zone to zone with the speed 
of light, and the kinetic one which has a 17-s rectangular 
shape delayed by 200 s and acts in the innermost zone of 
the shell. 

The shell mass exceeds _EshC~^ significantly (by a factor 
of 20 in our case). If the shell had uniform density distri- 
bution per unit solid angle, the shock would become non- 
relativistic well inside it, and there would be no 'canonical' 
long-term (several hours or days in observer's time-scale) 
power-law synchrotron afterglow. Thus, to produce such an 
afterglow, the shell should be inhomogeneous, consisting of 
'clumps' (or filaments) which intercept energy, and rarefied 
'windows', through which the relativistic fiuid could fiow 
without deceleration. As the shell itself is assumed to be cre- 
ated by a non-relativistic shock in the circumstellar medium, 
such 'clumps' and 'windows' may naturally form due to in- 
stabilities . Dense clumpy circum burst matter was also intro- 
duced bv lPostnov et al.l (|2004l ) to e xplain emission line s in 
the X-ray afterglow of GRB 011211 (JReeves et al.ll2002l ). 

Note, however, that for an opaque (or refiecting) inner 
boundary, during several days thermal emission can mimic 
the power-law decay Fi, oc t"'^ with exponential power typ- 
ical for observed GRB afterglows a ~ 1 — 2 (see Section 

4). 

Finally, some reduction procedures were added to take 
into account the jet beaming factor and the spherical shape 
of the shell: (1) the hotspot fraction is (1 — cos Sjot)/2 of the 
initially spherical shell; (2) the geometric time delay of the 
emission is 

5t{t) = (1 - e, ■ nob,)7?out(t - 5t{t))/c 

where Rout{t) is the radius of the outermost zone in which 
the thermal emission is calculated, e-r and nobs are local unit 
vectors directed from the emitting point radially and along 
the line of sight, respectively; (3) the angular distribution 
of the outcoming thermal emission (despite the main code 
is one-dimensional, an approximate angular distribution of 
radiation is repetitevily computed to update Eddington fac- 
tors). 

The maximal geometric curvature delay is 5t = (1 — 
Mjct)i?out(i — 5t)/c Ki 5000s, where fijet = cosSjet. The jet 
axis can be slightly misaligned with the line of sight, leading 
to a more complicated time averaging. Thus, for example, in 
Paper I a 3° off-set between these directions was assumed, 
resulting in the smearing off the curvature delay over 3000- 
8000 s interval, but in the present calculations no misalign- 
ment is assumed to see thermal features more clearly. 

On the illuminated spherical surface, a distant observer 
will see the emitting region as a set of concentric rings, each 
corresponding to the delayed surface luminosity. At a par- 
ticular moment, the observed isotropic luminosity equivalent 
will be the thermal luminosity in observer's direction aver- 
aged over the visible area with account of the time delay: 



ii',iso(t) = 87t^ / Iu{iJ.,t- 5t{t,^)')R^^^^{t ~ 5t{t,^)')^d^. 



First, a single-temperature fluid model is adopted for com- 
putational reasons. However, in a high-density medium (sev- 
eral 10^" cm~^), the full equilibration establishes in the flrst 
tens of seco nds (the co rresponding Uet parameter reaches 
10^^ cm~^s. IVinkI 120121 ). As the characteristic time-scales 
of the calculated light curves are significantly longer (thou- 
sands of seconds, see section 3), the one-temperature ap- 
proximation is justified. 

Second, we neglect multiple scatterings of gamma-ray 
photons. This allows us to recalculate the radiation field 
at each time step without considerable increasing of the 
computation time (see Paper I for discussion). Therefore, 
the total absorbed gamma-ray energy can be somewhat un- 
derestimated. For a given gamma-ray spectrum, the total 
heating power due to the second scattering is nearly four 
times smaller than due to the first one. Since the shell 
is semi-transparent for Thomson scattering, the account 
for multiple-scattered photons can not increase the heat- 
ing power significantly. However, it could affect outcoming 
gamma- and X-ray radiation spectra. 



3 RESULTS OF CALCULATIONS 

Fig. 1 presents some physical characteristics of the shell mat- 
ter during its heating: the gas temperature, its velocity and 
the thermal emission temperature Trad = (-4^^) 1 where 
Urad is the total energy density of the thermal emission. The 
profiles are shown for the time 800 s after the gamma-rays 
strike the shell. 

Fig. 2 demonstrates the light curves of the thermal emis- 
sion. Shown are the bolometric luminosity. X-ray (2-30 keV) , 
soft X-ray (0.1-2 keV) and UBVRI optical emission. For 
comparison, the characteristic range of the observed R-h&nA 
power-law afterglow light curves [j (corrected for the Galac- 
tic extinction and redshift) is also shown. Spectral density 
profiles for several characteristic moments are presented in 
Fig. 3 (e.g. the initial state, peaks of emission in different 
ranges; see the legend in Fig. 3). _R-magnitudes calculated 
for different redshifts (without extinction) are presented in 
Fig. 4. On these plots, time is counted from the moment 
at which the prompt emission front reaches the shell outer 
boundary, i.e. from the first pulse of the observed GRB. 

Gamma-ray effect. After the sharp front of the 
gamma-radiation with isotropic peak luminosity 3 x 10^^ erg 
s~^ impinges on the cold (3000 K) dense shell, the following 
sequence of events occurs. A nearly full ionization is reached 
in a few thousandths of a second since the photoionization 
rate dominates all other processes by several orders of mag- 
nitude. Only hydrogen-like ions of the heaviest elements are 
still present (0.1-10 ions per cm^, while the total particle 
number density is of order 10^^) due to recombination. This 
is important for the heat balance of matter and its thermal 
emissivity. 

As the material gets nearly fully ionized, the heating 
is mainly due to Compton scattering of gamma-rays off 
free electrons. The cooling is due to bremsstrahlung, inverse 



To calculate the light curves seen by the terrestrial ob- 
server, a cosmological model with Hq — 73.5 km s^^ Mpc~^, 
f^A = 0.76, r2m = 0.24 was used. 

Several simplifying assumptions should be mentioned. 



^ This range represents the upper and lower boundaries of the 
observed GRB afterglows, which we const ructed using about 4000 
data points c ompiled from the lit erature I IBadiin et al.ll2009r) . See 
also Fig. 7 in lGehrels et al.ll2009l for illustration. 
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Figure 1. Radial profiles of the fluid velocity V, temperature 
Tgas and the thermal emission temperature Tj-ad > calculated from 
its energy density, 800 s after the heating started. 



1======'^''''''''''''^'^^^'^'^'%^ 


^ 


— o — Bolometric 
—1— 0.1-2 keV 
^<— 2-30 keV 

U 

B 

V 
R 


\ \\ 


\ \ ^ 


' 





t 


V 



40 - 
39 - 



time, days 

Figure 2. Isotropic equivalent luminosities of the calculated 
thermal emission. Bolometric (50000A-30 keV, open circles), X- 
ray (2-30 keV, inclined crosses), soft X-ray (0.1-2 keV, upright 
crosses), UBVRI (solid, dashed, dotted, dash-dotted and dash- 
dot-dot lines, respectively) light curves are shown. For compar- 
ison, the range of the observed R-band power-law afterglows is 
plotted by the thin dotted straight lines. 



Compton scattering and radiative recombination which form 
the X-ray continuum. 

In 1 s after the radiative heating begins in a given 
fluid element, the dynamic balance of the medium and the 
gamma-radiation is reached at a temperature of ~ 10^ — 
10^ K, which varies as the gamma-ray luminosity changes, 
the sharp temperature spike in Fig. 1. 

After the main bulk of gamma-rays passes through a 
given fluid element, the partial recombination occurs pro- 
ducing mostly X-ray and soft X-ray photons. In a few min- 
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Figure 3. Spectral evolution of the thermal luminosity from the 
GRB-heated circumburst shell. The curves correspond to different 
stages (peaks, breaks etc., compare the given times with the light 
curves in Fig. 2) after the forward gamma-ray escape, i.e. after 
the GRB-trigger for an outer observer, except for the 'initial' one 
which shows the unilluminated shell spectrum. 



pushes it outwards with velocities of about 3000 km s~^ 
The outward thermal emission flux is accumulated near the 
gamma-ray front, while the inward flux maintains thermal 
balance in the underlying zones, somewhat decelerating their 
motion (the innermost zones may be even pushed back by 
this flux towards the centre of explosion), and then escapes 
through the inner boundary carrying the thermal energy out 
of the shell. 

GRB ejecta effect. After the jet kinetic energy has 
been injected (by means of the 'thermal bomb', see the pre- 
vious section), the medium gains an outward momentum 
and heat, and a radiative shock is formed, with the tem- 
perature initially rising up to 1.5 — 3 x 10® K (see the tem- 
perature bump of 10® K at R\xiS^i\ ~ 100.1 in Fig. 1, corre- 
sponding to the time around 600 s after the 'thermal bomb' 
had been activated). The dynamics of matter is essentially 
controlled by the radiation flux since the radiation pressure 
is several orders of magnitude as high as that of the gas 
even at the shock front. Thus, thermal radiation from the 
shocked matter pushes forward external zones before the 
shock reaches them, and the inward radiation flux (summed 
with that coming from overlying radiatively heated zones) 
decelerates underlying shell zones. Fig. 1 illustrates how the 
velocity jump is smoothed out and a cooling front moves 
from the transparent inner boundary into the shocked shell. 

As a result, the shock is smeared out and its strength 
is damped, so the outer zones coast with a velocity of 1.1 x 
lO'* km s~^ (in this model no external matter to the shell is 
assumed). 

Geometric curvature effect. While the local radi- 
ation field intensity changes on time-scales of ionization 



utes the temperature drops down to 8 — 10 x 10 K, and (~ 10 s), heating (several seconds) and cooling (minutes), 



stays approximately at this level until the gamma-ray front 
'causal discontiuity' crosses the shell and the thermal emis- 
sion begins to escape through the outer boundary (i.e. the 
temperature radial profile becomes rather flat behind the 
heating front). 

The momentum transferred by gamma-rays to matter 



the observed light curve will also be subjected to light travel 
effects: the shell light crossing time and the geometric cur- 
vature delay, which in our case are ~ 1700 s and ~ 5000 s, 
respectively. 

The shock emission becomes visible on the shell crossing 
time-scale (see optical curves in Fig. 2), while the bolometric 
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terrestrial observer's time, days 

Figure 4. Calculated ij-fluxes of thermal emission seen from dif- 
ferent redshifts. The straight dotted lines show the characteristic 
range of the observed _R-band afterglows of GRBs. 



luminosity maximum is achieved after about the curvature 
delay time, when the whole hot shell area is observed. After 
that, a characteristic time-scale of one or another light curve 
variability will reflect the cooling of matter in the shell (see 
e.g. a sharp break in the 2-30 keV curve in Fig. 2 indicating 
that the matter cools down to sub-keV temperature very 
quickly, during several seconds after the gamma-rays have 
escaped the shell and no longer supplied power to it). 

Thus, due to rapid gamma-ray heating and radiative 
cooling of the shell, a plateau with a steep decay emerges 
in the X-ray band. Thus, one should expect a more pro- 
nounced early and hard X-ray emission if the gamma-ray 
heating dominates. In the soft X-ray band a shallow (in log- 
arithmic scale) rise and shallow but stepening decay will 
be observed, which can be interpreted in terms of superpo- 
sition of the gamma-ray and the shock wave heating. The 
optical light curve is mainly shaped by the shock heating 
and thermal emission transport. This radiation should be 
more considerable if the most of the GRB energy comes in 
the jet kinetic energy. In extremal case where most of the ki- 
netic energy is thermalized in the shell, late-time bumps on 
the afterglow light curve may emerge, resembling a 'quasi- 
supernova' as discussed below. 

As can be seen in Fig. 4, the calculated R-fluxes of the 
optical thermal emission from sources at different redshifts 
lie within the region occupied by the observed R-afterglows 
of GRBs. The colour of the synchrotron optical spectrum 
is red, therefore the flux in the rest-frame R band of the 
source is higher than that in 7? band of the terrestrial ob- 
server. As the calculated thermal emission optical spectrum 
is bluer than the synchrotron one (see Fig. 3), for a terres- 
trial observer the thermal emission increases on top of the 
'red' synchrotron radiation, which facilitates its detection. 

The corresponding colour excess variations A(y — i?) = 
{V — R)tot — {V — R)phB are shown in Fig. 5 for different 
strengths and spectral slopes of the power-law background 
(denoted as PLB, Fi, ~ z^~ ) for the source at moderate {z = 
1) and high {z = 6.29) redshifts. The PLB spectral power 
indices are taken to be /3 = 0.7 with (V^ — -R)plb ~ 0.32 mag 
and /3 = 1.3 with (V - R)plb = 0.43 mag. The bright- 
est and the dimmest afterglow correspond to the upper and 
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Figure 5. Colour evolution of the thermal bump for different 
underlying afterglow power-law forms and source redshifts. A(V— 
R) is the change in the V — R colour relative to the PLB. The 
'bright' and 'dim' PLB stand for the upper and lower boundaries 
of the observed ij-afterglows in Fig. 4, respectively. The curves 
for bright afterglows are almost indisinguishable. 



lower boundaries of the observed power law afterglow range 
described above. 

As can be seen from Fig. 5, during the first hours the 
afterglow appears to be 'bluer' due to the contribution from 
thermal emission (although the total V — R colour is still 
positive), and its colour remains constant during the light 
curve bump. After the spectral maximum passes through 
the (blueshifted) R band, a noticeable reddening occurs. Ap- 
parently, the colour effect can be reliably observed only in 
afterglows with weak PLB. 



4 COMPARISON WITH OBSERVATIONS 

Optical range. Optical afterglows of several GRBs demon- 
strate clear deviations from the pure power-law decay 
(see e.g. discussion in Paper I). These 'bumps' correspond 
to additional proper luminosities of the order of 10'*^ — 
10*^ erg s~^, which is close to what is expected from thermal 
emission of circumburst shells illuminated by prompt GRB 
emission. 

Because of cosmological 7^-corrections the luminosity 
spectral density appears to be not a good estimator to com- 
pare observations and calculations due to cosmological K- 
corrections. Instead, we use the observed 7?-band magni- 
tudes, which can be derived from published data (corrected 
for the Galactic extinction), and compare them with model 
i?-magnitudes, which can be readily obtained as {1 + z) 
times shifted flux spectral density convolved with the R- 
filter transmission function. 

Table 1 lists properties of several GRBs which we use 
for comparison with our model. It shows that the observed 
bumps occur close to the time expected for thermal emission 
from the illuminated shell. For GRB 050904 no spectral in- 
dex is available (it was not observed in R-Hltei because of its 
high redshift and the host galaxy Lyman-break absorption 
falling into i?-band). The spectral slope and extrapolated 
7?-magnitudes (ignoring the host gal axy absorption) are es- 
timated using available NIR colours (|Tagliaferri et al.ll2005l : 
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Table 1. Paramete rs of bumps of optical afterglow of GRBs obtained from data collected in 
iBadiin et al.1 1 120091 ) and corrected for the Galactic extinction. 
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Figure 6. The GRB 020124 optical afterglow in the R band (open 
circles) and the 2 = 3 calculated thermal emission i?-band stellar 
magnitudes (solid line). 
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Figure 7. GRB 050904 optical afterglow extrapolated to R band 
(open circles) and the z = 6.29 calculated thermal emission stellar 
magnitudes in the R (dashed line) and / (solid line) bands. 



IPrice et al.1 12 006: Hais lip et all |2Q06| : iBoer et al.1 12006| . see 
also lKann et al. 2006). 

Fig. 6 and 7 demonstrate 7?-light curves of GRBs 020124 
and 050904 (synthesized) with the calculated thermal emis- 
sion superimposed. Clearly, the agreement is not perfect (we 
have not attempted to tune our model to fit any specific fea- 
ture), but the similar time of emergence and amplitude of 
thermal emission from the gamma-ray illuminated circum- 
burst shell in both cases seems to be suggestive. 

As can be seen from Fig. 6 and 7, the thermal emission 
does not cor rectly reproduce t he on set of the irregularity. 
However, as iNakar fc Granoj l|2007l ) show, such a smooth 
transition without a considerable rebrightening could be due 
to contribution of the shock synchrotron radiation when the 
shock encounters a sharp density contrast. The transition 
should become visible at times of several To (the time when 



the shock strikes the dense region), i.e. on the time-scale 
of several (5f^_ah multiplied by a relevant (1 -I- z) factor. 
The adopted 5t^,-sh ~ 10 s, thus, gives appropriate times 
(since GRB trigger) of the irregularity onsets of ~ 0.1 d. 
Later, when the thermal emission gains strength and dom- 
inates, it forms long-term bump and rapid fading (which 
|Nakar fc G ranot 200'^ argue to be not produced by the syn- 
chrotron emission only). 

Note also that all considered GRBs have relatively low 
spectral index (3, i.e. they tend to be systematically 'bluer' 
than the typical power-law afterglows with mean /? ~ 0. o 



■^ The spectral slopes have not been separately derived for the 
bumps; th ey have been taken from observational reports (see ref- 
erences in lBadiin et al.ll2009h . and hence are averaged over time. 
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Figure 8. The GRB 050904 (grey crosses) and GRB 070110 
(black crosses) X-ray luminosities obtained bv lTroia et alj BOOT! ) 
from observed 0.3-10 keV fluxes, and model thermal emission light 
curves: bolometric (open circles), 0.1-2 keV (dashed line) and 2-30 
keV (solid line). 



This also does not contradict their interpretation in terms 
of thermal emission from circumburst shells. 

X-rays. The calculated spectral density of thermal 
emission (see Fig. 2) suggests that most of the shell ther- 
mal energy is carried away by X-rays. The properties of the 
emergent X-ray emission are directly sensitive to parameters 
of the shell and heating conditions, while the dependence of 
the optical emission on the shell parameters is less straight- 
forward to explain. The most notable feature of the thermal 
emission in X-rays is the appearance of a 'plateau' or slowly 
rising bumps mainly due to the high-latitude photon time 
dela y from the shell's illuminated part. 

Troia et al.ll2007l present observed 0.3-10 keV afterglows 



of GRB 050904 and 070110 converted into luminosity ver- 
sus proper time, demonstrating a plateau exactly where the 
thermal bump is expected. In Fig. 8 we have placed both 
these light curves and our calculated bolometric. X-ray and 
soft X-ray ones. 

A thermal bump from the illuminated shell provides suf- 
ficient power to explain the feature observed in GRB 070110 
afterglow, but underestimates the amplitude of the X-ray 
plateau observed in GRB 050904 by an order of magni- 
tude. However, GRB 050904 itself was intrinsically brighter 
by an or der of magnitude according to different estimates 
(see e.g. iKann et al.l 120061 and references therein). The X- 
ray plateau in GRB 070110 had a flux of about 2 (i.Jy at 
1 keV corresponding to a luminosity spectral density of 
Ly ~ 5 X 10^^ erg s~^ Hz~^ at 3.3 keV in the source frame, 
which is in good agreement with thermal emission model 
predictions (see Fig. 3). The calculated light curves also re- 
produce the abrupt end of both the plateaus occurring at 
the same time. 

A certain discrepancy in t he hardness ratio (1-10 keV 
to 0.2-1 keV) should be noted. iTroia et all IJ2007h reported 
that the observed hardness ratio of GRB 070110 was about 
one and increased up to ~ 1.5 during the plateau phase, 
while the main part of our model thermal emission comes in 
soft X-rays (0.1-2 keV). 

The deficit of hard photons in our model can be par- 



tially due to the one-temperature approximation. Indeed, 
the nearly instantaneous photoionization injects a lot of elec- 
trons with energies up to several ionization energies of ions 
they were knocked off, i.e. up to tens of keVs (~ 10* K), 
much higher than the ion temperature prior to the gamma- 
ray flashing. Compton scattering can additionally heat them 
up to 60 keV (the Compton temperature corresponding 
to the adopted gamma-ray spectrum and electron density, 
though this is only an approximate estimate, since the heat- 
ing and cooling are non-stationary). Thus there should be 
more bremsstrahlung in X-ray and hard X-ray bands before 
they cool down to the equilibrium with ions (lO'' K). 

Compton-scattered prompt emission photons should 
also noticeably contribute to hard X-rays and increase the 
hardness ratio. The prompt emission 1-60 keV photons scat- 
tered towards the observer can provide an X-ray luminosity 
of L^60,iso ~ 1.25 X 10"*^ erg s~^ (assuming the bolometric 
luminosity 3 x 10^^ erg/s). Also, some part of the inward 
coming thermal photons should be upscattered into hard X- 
rays. The scattered emission will have a non-thermal spec- 
trum harder than that shown in Fig. 3, and its light curve 
will show a rise on the curvature delay time-scale (i.e. thou- 
sands of seconds) followed by abrupt fading. This seems to 
be very similar to what is needed to produce the hard X- 
ray plateaus in the light curves of GRB 050904 and GRB 
070110. 

Thus, the observed hard X-ray irregularity, which can 
have a non-thermal origin, may indicate a shell-like circum- 
stellar feature, which in turn may produce a thermal emis- 
sion detectable as another bump in softer band at later time. 
Therefore we argue that the similarity between the observed 
afterglow bumps and the calculated thermal emission prop- 
erties (in magnitudes, characteristic times, 'bluer' spectral 
slopes. X-ray features) may be a strong indication of the 
presence of dense circumburst shells around some GRBs. 
This should be taken into account in studies of GRB after- 
glows. 

Effect of boundary conditions: the 'quasi- 
supernova' effect. When the inner boundary of the shell 
is rather opaque (or reflective, if there is a back scattering at 
the relativistic shock), the radiation does not escape through 
it from the shell at all or diffuses slowly, while it can freely 
escape outwards. During the radiative heating phase (when 
the intermediate causal discontinuity exists) the energy is 
trapped within hot parts of the shell (between the prompt 
emission front and the inner boundary) resulting in strong 
radiative forces. It dominates also after the energy injection 
('thermal bomb') due to the impact of the ejecta. During 
the flrst hours all motions of the hot medium are controlled 
by thermal radiation fluxes directed outwards. 

The main differences with the model considered above, 
in which the inner boundary was transparent, are as fol- 
lows. The 2-30 keV luminosity and its characteristic time- 
scale have not changed significantly, but in soft X-rays and 
UBVRI bands a bright flash appears produced by the shock 
impact, which becomes nearly ten times more powerful than 
the radiation produced in the case with the transparent in- 
ner boundary. 

This is followed by a several-day-long fading with the 
slope similar to the ordinary power-law afterglow. During 
this interval, there is no or only weak optical colour evolu- 
tion, therefore the sum of thermal and non-thermal emission 
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Figure 9. Calculated ij-fluxes of the 'quasi-supcrnova' thermal 
emission from different redshifts. The straight dotted lines show 
the observed i?-band afterglow range. 



will behave as an ordinary afterglow (possibly slightly bluer 
than the usual one). 

The shock and the radiation pressure accelerate the 
shell up to 5 — 6.5 x 10* km s~^, so several days after 
the explosion, a bright supernova-like bump emerges, with 
colour changing from blue to red in a way similar to what 
is observed in a ty pe Iln supernova originating in a dense 
circumstellar wind l|Filippenkdll997l ). Model optical bumps 
observed from different redshifts are shown in Fig. 9. 

Interestingly, similar double-bumped light curves were 
observed by gw/i-UVQT in GRB 060218 (SN2006aj) af- 
terglow by ICampana et al.l (|2006r ). A soft X-ray ther- 
mal component slowly rising for 3000 s and then fading 
abrubtly was also detected (see Fig. 2 in their paper). 
Rejecting a criticism cast on the feasibility of the shock 
breakout to produce a bright enoug h thermal emission, 
IWaxman. Meszaros &: Campanal (|2007l ) suggested the pres- 
ence of an additional circumstellar structure (envelope) in 
the dense wind. 

Another evidence for a dense G RB environment was 
discussed by (|Garnavich et al.l 120031 ) . who noted that the 
SN2QQlke assigned to the GRB 011121 was of type Iln rather 
than Ib/c. 

Therefore, the radiation-driven 'quasi-supernovae' 
might be responsible for late-time bumps in optical after- 
glows of GRBs, which are usually classified as supernova- 
bumps. However in the present calculations, it appears 
about 4-5 magnitudes brighter than required to explain 
these bumps, and more detailed treatment of the inner 
boundary condition as well as non-spherical symmetric ra- 
diation transport is needed. 

Either a strong opasity contrast or refiection is required 
to produce high radiation flux gradients driving the 'quasi- 
supernova'. For example, high radiation flux gradients can 
appear due to high density contrasts in regions of the ex- 
plosion, implosion, shocks or cumulative jet developement, 
where the radiation or the ejecta propagate in a conical 
channel encircled by a cold dense medium, or there are 
clumps of matter on the way, when the jet protrudes through 
the progenitor star or its extended atmosphere, when the 



blown-away material expands into the stellar wind, or when 
the GRB radiation or ejecta perturb a filament of the cir- 
cumburst matter, etc. Additionally, as mentioned above, 
the relativistic aberration of the emission scattered within 
the ejecta can effectively act like reflection. In such cases 
supernova-like features can be observed in the afterglow, 
without supernova produced by the exploding star itself. 

It is tempting to apply this model to the explanation 
of the GRB-SN connection. While the GRB generating sce- 
nario is fr equently referr ed to as a coUapsar or a 'failed su- 
pernova' (jWooslevI Il993r ) . a fairly 'successful' supernova is 
required to produce observed late-time afterglow bumps. If 
both phenomena originate in the progenitor star interiors, 
then the same GRB central engine must be able to drive 
both the narrow ultrarelativistic outflow with a low bary- 
onic load, and the quasi-spherical explosion with a kinetic 
energy of the same order. Moreover, the GRB supernovae 
are found to be systematically intrinsic ally brighter tha n 
their regular mates (see e.g. section 4 in lCano et al.l 120111 ). 
However, this difficulty might be overcome if the central en- 
gine launches the jet only, and the source of the observed 
supernova 'bump' lies outside the central star. In this case, 
the GRB supernova would be induced by the deposition of 
a portion of the GRB kinetic and/or radiation energy into 
the progenitor star exteriors or its environment. 

Clearly, the problem of the GRB jet/shell interaction 
and the emerging radiation is multidimensional and rela- 
tivistic, and so far no self-consistent solution has been ob- 
tained. Researchers working in the field either simplify (or 
neglect at all) the radiation transport, d eriving the radiation 
properties from the h ydrodynamics fe.g. lLazzati et al.ll2009l : 
iNagakura et al.ll2011^ ■ or limit themselves only by a simpli- 
fied hydrod ynamics (like the Blandford-McKee self- similar 
solution, see lTolstov fc Blinniko-^120031 : lTbrstovll2010l ). 

We emphasize that thermal emission from circumburst 
medium can significantly contribute to the radiation of 
GRB afterglows. The produced thermal radiation is strongly 
coupled with properties of non-thermal emission of GRBs. 
Therefore, the development of sophisticated multifrequency 
multidimensional relativistic radiation hydrocodes to treat 
interaction of radiation with matter is strongly desirable to 
improve our understanding of the GRB physics. 



5 CONCLUSIONS 

In the present paper we have used the modified radiation hy- 
drocode STELLA to calculate thermal radiation effects from 
the interaction of the prompt GRB emission and GRB jet 
with the dense shells of matter which can surround mas- 
sive GRB progenitors. Our model calculations revealed that 
thermal radiation of the heated circumburst structures can 
be visible in the X-ray and optical GRB afterglows in the 
form of plateaus and deviations from the pure power-law 
time decay. In particular, the calculated thermal radiatiton 
from a dense shell with the total mass of several Mq lo- 
cated at a distance of 10^^ cm from the GRB centre, which 
is illuminated by a 17-s long GRB with the total energy 
4.5 X 10^^ ergs, can reproduce both the time, shape and 
magnitude of real X-ray plateaus and irregularities of opti- 
cal afterglows in several GRBs. It is found that the coupling 
between the radiation and dynamics of matter in such shells 
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can lead to an interesting phenomenon - a 'quasi-supernova' 
effect (radiation-driven explosion of the shell). This intrigu- 
ing possibility can explain the GRB-SN connection without 
requiring in some cases the quasi-spherical explosion of the 
progenitor star itself, and is worth further studying by means 
of more realistic multi-dimensional relativistic radiation hy- 
drocodes. 
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